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Abstract 

We obtain a limit when mass tends to zero of the relativistic difFusion of 
Schay and Dudley. The diffusion process has the log-normal distribution. 
We discuss Langevin stochastic differential equations leading to an equi- 
librium distribution. We show that for the Jiittner equilibrium distribu- 
tion the relativistic diffusion is a linear approximation to the Kompaneetz 
equation describing a photon diffusion in an electron gas. The stochastic 
equation corresponding to the Jiittner distribution is explicitly soluble. 
We relate the relativistic diffusion to imaginary time quantum mechanics. 
Some astrophysical applications (including the Sunyaev-Zeldovich effect) 
are briefly discussed. 

1 Introduction 

An extension of the idea of difFusion to relativistic theories poses some problems 
of conceptual as well as physical character (for a discussion of these problems and 
different approaches to solve them see ([1] [2] [3] [4] [5] [6]; for a review and further 
references see [7] [8]). In [9] we developed the difFusion theory of Schay [1] and 
Dudley [2] of massive relativistic particles discussing friction terms leading to an 
equilibrium. In refs.[l][2] the relativistic Brownian motion is uniquely defined 
by the requirement that this is the diffusion whose four momentum stays on the 
mass-shell. The formulas become singular as m — > 0. 

The photon is a massless particle which could experience a diffusive behav- 
ior. We could think of a multiple scattering of the photon in a charged gas as 
resembling the Brownian motion. The number of such encounters must be large 
if the approximation is to be feasible. This happens in astrophysical applica- 
tions when the photon travels over large distances in intergalactic space filled 
by a ionized gas. The diffusive approximation could also apply to the early 
stages of the Big Bang and to the photon motion inside the plasma shell of the 
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compact stellar objects (stars, black holes) The relativistic Boltzmann equa- 
tion is the standard tool in a study of the evolution of the photon in a gas of 
charged particles [10] [11]. The diffusion approximation to the Boltzmann equa- 
tion has been applied in many models of the transport phenomena [12]. In the 
context of the photon propagation such a diffusion approximation is known as 
the (non-linear) Kompancctz equation [13] [14] [10]. The Kompancetz equation 
arises as an approximation to the Boltzmann equation describing an evolution 
of the photon in a gas of electrons. In such a case the evolution of a beam of 
photons is disturbed by Compton scattering and Bremsstrahlung. These are 
the processes widely studied in astrophysics [10] [15] [16]. 

In this paper we investigate a diffusion of massless particles without spin 
in the framework of classical relativistic diffusion theory . We show that there 
is a massless counterpart of Schay [1] and Dudley [2] relativistic diffusion (the 
formula for a diffusion of a massless particle has been also derived in [17] as a 
continuum limit of discrete dynamics of space-time, sec also [18]). We derive 
the diffusion as a classical limit of quantum mechanics (the Boltzmann equation 
usually also incorporates quantum scattering processes). Such a diffusion pro- 
cess has no limit as time grows to infinity. We discuss (sec. 3) drags which lead 
to an equilibration of the process. Subsequently, we show in sec. 4 that if the 
Kompaneetz equation is expanded around its equilibrium solution (Jiittner or 
Bose-Einstein) and non-linear correction terms are neglected then the resulting 
diffusion equation is the same as our relativistic diffusion with the friction lead- 
ing to the same equilibrium. In sees. 5 and 6 we discuss the dynamics in an afRne 
time parameter. It is shown that the dynamics (in particular, the approach to 
the equilibrium) can be well controlled by methods of Langevin theory (sec. 5) 
or quantum mechanics (sec. 6). In the final sec. 7 we discuss the problem of time 
in the relativistic diffusion. We show that it can be approached in a similar 
way as in the deterministic relativistic dynamics. The af&ne parameter on the 
trajectory can be replaced by the (random) laboratory time. As a result from a 
solution of the diffusion equation in the afhnc time we can obtain a solution of 
the transport equation involving only the physical variables (in particular, the 
laboratory time). 

2 Quantum origin of the relativistic diffusion 
equation 

In our earlier paper [9], following Schay [1] and Dudley[2], we have defined the 
generator of the relativistic diffusion as the second order differential operator 
on the mass-shell Hm 

= Voirf - pi{Tf -p2iT)^ -P3(t)^ ^m^c" (1) 

The generator is singular at m = 0. For this reason in order to define a diffusion 
of massless particles we propose another method. A diffusion should be con- 
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sidered as an approximation of a complex multi-particle dynamics. In our case 
the relativistic dynamics. It is not clear how to describe many particle systems 
in a relativistic way. Quantum field theory leads to a relativistic description of 
scattering processes. However, its reduction to particle dynamics has not been 
explored. The photon as a massless particle is a typically relativistic object. At 
the same time this is a quantum system with an internal angular momentum 
(spin). We should describe it by a relativistic wave function. In this paper we 
neglect the spin. In [19] we discuss the diffusion of particles with a spin (the 
helicity in the massless case). It is shown in [19] that the dissipative part of the 
evolution does not depend on the helicity if m = 0. Hence, the neglect of spin 
in this paper is justified. 

The explicitly Lorentz invariant description of dynamics must in fact be 
static. It has to be described in terms of space-time trajectories. The laboratory 
time a;*^ is one of the coordinates. Hence, the evolution as a function of this 
coordinate should be inferred at the later stage from a study of the set of paths. 
In the massive case the position of a particle on the space-time path is most 
conveniently described by its own time : the proper time. If a particle is massless 
then the notion of the proper time as a time of an observer moving with the 
particle does not make sense. Nevertheless, it is still convenient to introduce an 
affine parameter r describing the fictitious motion along the particle trajectory 
(we shall still call it a proper time although no observer can move together with 
the particle). 

The proper time as an additional parameter in a description of a relativistic 
quantum particle has been introduced in [20] [21] (in quantum theory the relation 
to the time of an observer moving with the particle is lost anyhow). It can be 
treated as a convenient tool in a formulation of relativistic quantum mechanics 
and quantum field theory. The equation for a massless free field is 



where 9q — A is the wave operator. 

Wc may consider a perturbation of (2) by some other fields (an environment) 
so that eq. (2) remains invariant under the Lorentz group. When we average over 
the environment then in the Markovian approximation the system will be de- 
scribed by the master equation for the density matrix p. A preservation of the 
trace and positivity of the density matrix requires that the master equation 
should have the Lindblad form [22]. Then, the Lorentz invariance and an as- 
sumption that the Lindblad generators arc built from the Lorentz generators 
leads in the simplest (linear) case to the equation 



where are the generators of the algebra of the Lorentz group. Eq.(3) can 
be rewritten as an equation for the Wigner function [23] (proper time equations 



id^cj) = {dl - A)(j) 



(2) 



idrP = [dl -A,p] + \i^[M^., [M''^ p]] 



(3) 
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for the Wigner function in relativistic quantum mechanics appear in [24] and in 
quantum field theory in [25] [26]) 

where 

V = -^(P.^-f^^) (5) 

is a reahzation of the algebra of the Lorentz group in the momentum space. 
It can be checked that the diffusion generator of [9] can be expressed as 

Ah = ^i,.^^''" (6) 

The derivation of the relativistic diffusion equation from quantum mechanics is 
similar to a derivation of the Boltzmann equation from quantum mechanics in 
[27] . Eq. (4) makes sense in the massless case. Choosing p as coordinates on Ho 
we obtain 

Ah = PjPkd^d'' + Sp'^dk (7) 
where k = 1,2,3 and = We note that the generator (7) is the limit 

m ^ of m?c^/\H of rcf.[9]. 

The relativistic diffusion equation in the momentum space (a relativistic 
analog of the Brownian motion) reads 

dr(l)r = (8) 

7^ is a diffusion constant which has the dimension of (our choice in the next 
sectionrmomentum divided by length). 



3 An approach to the equihbrium 

We add a drag term 



to the diffusion (8). The vector field (9) defines a perturbation of the relativistic 
dynamical system (at 7 = 0) 
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where from eq.(l) (for m = 0) we have p° = |p|. The zero component of 
eq.(ll) determines a relation between the parameter r on the trajectory and 
the laboratory time . Let 

g='^^''AH + X (12) 
A solution of the diffusion equation 

dr(l>r = Q(t>r (13) 

with the initial condition 0(x, p) defines a distribution of random paths (see 
sec. 5) starting from the point {x, p) of the phase space. At 7 = we have 
a distribution of deterministic trajectories. The dynamics in the laboratory 
time must be derived from the static picture of space-time paths ( see the 
discussion of relativistic statistical dynamics in [3]). 

We are interested in the behavior of the relativistic dynamics at large values 
of T. In particular, whether the diffusion generated by Q can have an equilibrium 
limit. We say that the probability distribution dxdp$ is the invariant measure 
for the diffusion process (see [28]) if 

j dxdp^{x,p)(pr{p,x) = J dxdp^r{p,x)^{p,x) = const (14) 

where 



dr^r = g*<^r (15) 

and g* is the adjoint of Q in L'^{dpdx) (the Lebesgue measure dp is not Lorentz 
invariant, the factor |p|~^ necessary for the Lorentz invariance is contained in 
see [9]). Hence, the invariant measure dxdp^R is the solution of the transport 
equation g*^R = 0. Exphcitly, 

^^^dW%pk^R - ^j^dipj^R - ^^dmj^R = p^dli^R (16) 

(derivatives over the space-time coordinates will have an index x). We denote 
by ^ER the x-independent solution of eq.(16). Then, R can be expressed in 
terms of ^er 

Rj = -^P] + ■^PjVkd'' In ^ER (17) 

We assume that ^er is a function of the energy pqc multiplied by a constant 
(3 of the dimension inverse to the dimension of the energy (/? = in the 
conventional notation). In such a case from eq.(17) 

Ri = \p, + \pcp,\p\{\n<^ER)'mp\) (18) 
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Now, eq.(13) reads 

dr<j>r = h^Ancl^r + ll'vA^ + /3c|p| (In (19) 

and eq.(16) 

l^-^did^VjVk'^R - 2^^d^Pj^R - ^j^pcd^pj\p\{ln ^erY^r = P^^dl^^R (20) 

Let us consider an initial probability distribution of the form $ = ^er^- Then, 
from eq.(15) we obtain its evolution 

= ^ER^r (21) 

as a solution of the equation (we do not make any assumptions on the form of 
^ER here ) 

(22) 

4 Spherical coordinates and the Kompaneetz equa- 
tion 

In this section we restrict ourselves to the diffusion in momentum space. The 
diffusion generator(7) is degenerate. This is easy to see if we express it in the 
spherical coordinates on Ho 

Po = r (23) 
Pi = rcos(psm6, P2 = rsm^sm6,p3 = rcos9. In these coordinates 

AH = r2|;+3r# = #; + 2# (24) 

where we introduced an exponential parametrization of r 

r = exp u (25) 

Here, u varies over the whole real axis. It follows that the three-dimensional 
diffusion of massive particles becomes one-dimensional in the limit m ^ 0. 
In the coordinates (23) it is sufficient if we restrict ourselves to the drags 

X = ,^rI (26) 
The diffusion (12) is generated by 



g = ^idl + 2d^)+j^Rdu (27) 
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We may write for integrals of spherically symmetric functions 

dp = 4:TTdu exp(3u) (28) 

Then 

g* = e-^"g+e^" (29) 

where 

g+ = ^{dl-2du)-'y^d^R (30) 

is the adjoint of Q in L^{dudx). The probability distribution evolves according 
to eq.(15). The invariant measure dxdp^u = dxdu^j solves an analog of eq.(16) 

g*^R = g+^i = (31) 

where 

$/ = exp(3w)$i? (32) 
We can express eq.(31) as an evolution equation in the laboratory time 

5o$/ = exp(-?i)pVx$/ 

+^exp(-u)(a2-2a„)$/-72cxp(-«)a,(7?$7) ^ ' 

where n = exp(— u)p is a unit vector in the direction of motion (equations of this 
type describe a diffusive perturbation of the propagation of a massless particle 
with the velocity of light [29]). 

The independent solution of eq.(33) is denoted The drift R is related 

to 

R=^{du\n^Ei-'i) (34) 

If the diffusion (15) starts from an initial distribution 

$/ = $B/*^ (35) 

then from eq.(15), using the equilibrium relation (34), we obtain the evolution 
of* 

dr^i = li'^dl^l + \^\du \n<^Ei)du^i (36) 
We consider Jiittner equilibrium distribution [30] in L^(du) 

^Ei = exp(-/?cr) = r^nj (37) 

Then, from eq.(34) 

R = -^{-1+ (3cexp{u)) (38) 
The diffusion generator corresponding to the Jiittner distribution reads 
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For the Bose-Einstein distribution (in Compton scattering the photon number 
is preserved; if an equilibrium is achieved, then the chemical potential ^ 0, 
see [14]) 

$1^ = ( exp(/3(M + cr)) - l) = r^UEif^) (40) 

we have 

-R = -^7^1^ - 1 + /?cexp(w)(^l -exp(-/3/i-/3cexpu)^ ^ (41) 

Eq.(36) for an evolution of '^^ (with the drift (41) leading to the Bose-Einstein 
equilibrium distribution) in the r-coordinates reads 

9^4-^ = ^j^r-^(dydr'^^ + r\dr\nnE)dr'^^) (42) 

where 

drlnnE = -(3c{l - cxp(-/3^ - /Scr))"^ (43) 

In the low temperature limit (or for the Jiittner distribution with the generator 
(39)) we have in the diffusion equation (42) 

dr InnE — > —/3c 

We compare the relativistic diffusion equation with the (non-linear) Kom- 
paneetz equation usually written in the form [13] 

drn = K^x-'^d^(x'^{d^n + n + n^)j (44) 

Eqs.(42) and (44) coincide if the non-gradient and friction terms arc neglected. 
Note that the form of the differential operators on the rhs of eqs.(42) and (44) 
follows from the conservation of probability for eq.(42) and conservation of the 
number of photons for the Kompaneetz equation, i.e.. 



In order to explore the appearance of friction in cq.(44) we expand the photon 
distribution of the Kompaneetz equation around its equilibrium value in the 
same way as we expanded the distribution of the diffusion process (15) around 
the equilibrium in eq.(36). If we neglect on the rhs of cq.(44) then the rhs 
disappears (because dxnj = —nj, with x = Per) for the Jiittner distribution. 
Let 

n = exp{—x)'if 

Then, neglecting the n? term in the Kompaneetz equation (44) we obtain (after 
an elementary rescaling r f3cr = x, = ^7^) the diffusion equation (42) 
with Ue Tij. 
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In general, the Kompaneetz distribution n equilibrates to ns (as dxUs + 
riE + n% = on the rhs of eq.(44)). Let us write the initial distribution in the 
form 

n = heT 

then the evolution of F is determined by the equation 

(45) 

The last term in eq.(45) is absent in eq.(42). It describes the effect of the 
quantum statistics which promotes bosons to condense. As long as the density 
of photons is low this term will be negligible. 

It remains to consider the Lorentz transformations of the equilibrating diffu- 
sion equation. The diffusion equation (8) (without friction) is explicitly Lorentz 
invariant as follows from eq.(6). However, the equilibrium measure cannot be 
Lorentz invariant if it is to be normalizable (in a finite volume) and if the mass 
shell condition (1) is to be satisfied (there is only one invariant measure on the 
mass-shell, this is dppQ^). The equilibration to the Jiittner or Bose-Einstein 
equilibrium takes place in a preferred Lorentz frame [31]. In a general frame 
described by the four-velocity w" we should write the equilibrium distribution 
in the form ^ 

nE{w,fj,) = (^ejip^Pw^p^ + Pw^^ity) - 1^ (46) 

where (by convention) in the rest frame w = (1, 0, 0, 0) and /xq = /i. 
We can write the Kompaneetz equation (44) in an arbitrary frame 

drn{w) = ^(h^pL'^Pniw) 

\ (47) 

+2iw^L''Pppn{w){n{w) + 1) + 2w''p^n{w){n{w) + 1) j 

In the rest frame eq.(47) coincides with eq.(44) as can be checked using eq.(5) 
(with Loj = —i\p\dj on the mass shell). In an arbitrary frame nr{w) has the 
Bose-Einstein distribution nEiw,iJ,) as an equilibrium measure. 

The Kompaneetz equation finds applications in astrophysics [10] [32] [33] 
[34] [35] [36]. The CMBR photons described by n£;(/i = 0) are scattered 
when passing clusters of galaxies and other reservoir of hot plasma. Eq.(45) 
could be solved with ueIp- = 0) and F = 1 as an initial condition (see the 
end of the next section). Then, it describes the photon density evolution as 
it passes through the electron gas. The resulting distortion of the black-body 
spectrum (Sunyaev-Zeldovich effect, see [34] [37] for reviews)) can be compared 
with observations (as discussed in [19] taking into account the photon helicity 
does not change the diffusion equations (8) and (13), the spin dependent part 
disappears in the limit m ^ 0). The Kompaneetz equation has been derived 
from the Boltzmann equation [13] [10] under an assumption that the electron ve- 
locities are non-relativistic. Relativistic corrections to the Boltzmann equation 
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describing the Compton scattering in a gas of electrons and its diffusion limit 
have been calculated in [16] [35] [38]. These computations involve relativistic cor- 
rections to the photon-cIcctron cross section in the laboratory rest frame. They 
are unrelated to the generalization (47) of the diffusion equation to an arbitrary 
frame. 



5 Langevin stochastic equations 

We wish to express the solution of the diffusion equation (8) by a stochastic 
process 

4>r{p,x) = E[(^{p{t),x{t))] (48) 

Here, 



x^'{T) = a;" + / p^{s)ds (49) 
Jo 

The diffusion process p{t) as well as its generator Q can be expressed in various 
coordinates. In the coordinates (7) we have to find a frame Cja {j,a = 1,2,3) 
such that 

PjPk = GjaBka (50) 

A solution of eq.(50) is 

Sja = Pjda (51) 

where CaCa = 1- Then, the stochastic process solving eqs.(8) and (48) reads 
[28] [39] 

3 

^Pj = + IPj^adh" = ^"^Pjdr + ^PjSa o db" (52) 

where b" are independent Brownian motions defined as the Gaussian process 
with the covariance 

E[b''{T)b''{T')] = S'^'miniT, t') (53) 

(we denote the expectation values by £^[...]). The first of eqs.(52) is the Ito 
equation whereas the circle in the second equation denotes the Stratonovitch 
differential (the notation is the same as in [28]). Eq.(52) is equivalent to the 
simpler equation for r = jpj 

dr = j'^rdr + 'yrodb (54) 
where b is the one-dimensional Brownian motion. The solution of eq.(52) is 

Pji'^) = Pj exp(7^T + 7ea6"(T)) (55) 
whereas the one of eq.(54) 

Tr = r exp(7^r + 76(t)) (56) 
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We could derive eq.(56) from eq.(55) setting b = Cab""- Eq.(56) means that (in 
the coordinates (25)) 

Ur = u + 7^r + 76(t) (57) 

When we know the probabiHty distribution of then we can calculate the 
probability distribution of the momenta. We have (in the coordinates (23) the 
angles do not change in time) 

E[f iPr)] = E[f{rr, <P, 9)] = (27r7V)-i /~ f /(y, </., 9) exp ( - ^(In ^ - j^r 

Hence, IpI has the log-normal distribution. Note that if we wish to express 
the probability distribution "if (42) when drlnns in terms of a stochastic 
process Ur then Uj- is the solution of the stochastic equation 

3 

dur = 'jl'^dr + 'ydbr 

In such a case 

*.(r) = i?[*(exp(«.))] = (27r7V)-i |^4'(2/)exp ( - ^^(Inf - §7^)^) 

(58) 

is the solution of the equation 

with the initial condition ^. The solution (58) is discussed in [33] [36]. For small 
T and the initial condition ^ = nsilJ' = the solution can be approximated by 

*r = nB(M = 0) + Tg*nE{n = 0) 

This approximate solution of the diffusion equation without friction is usu- 
ally applied in a discussion of Sunyaev-Zeldovich CMBR spectrum distortion 
[33] [34] [35] [36] [37]. 

The stochastic equation for the diffusion with the Jiittner friction (38) gen- 
erated by (39) reads 

du = j^dr — ^7^/?cexp(u)(ir -|- jdbr (59) 

It has the solution 

Ur = u + j'^T + jbr - In ^1 + i7^/3cexp(u) JJ" exp{'ybs + 7^s)rfs^ (60) 

It is not easy to calculate the correlation functions of Ur or in a closed 
form (we could do it as an expansion in /3). However, correlation functions of 
^ = exp(— u^) can be calculated explicitly. 
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6 A relation to the imaginary time quantum me- 
chanics 

We can express solutions of the diffusion equations by an imaginary time evolu- 
tion generated by a quantum mechanical Hamiltonian.The diffusion generator 
is of the form 

Q=\l'Ql-^Qu (61) 



2 



Let 



Then 



where 



(62) 



exp{-n)g exp{n) = -H= ^^"^dl - V (63) 



We have the standard Feynman-Kac formula [28] [39] for exp(— rif) 
(exp(— Ti?)V')(w) = E exp - j V{u + 'ybs)dsjtjj{u + 76^) 



(65) 



Let 

ip = exp{-Q)(t> (66) 

Then, from eq.(63) it follows 

4>riu) = iexp{Tg)(t>){u) = E[<f>{ur{u))] = exp(0(u))(exp(-Tfl')V)(w) (67) 

where Ur (u) (for the Jiittner distribution) is the solution (60) of the stochastic 
equation (59) with the initial condition u. Inserting ^ = 1 in eq.(67) we obtain 
exp(— ri?) exp(— fi) = exp(— fi). Hence, 

Hexp{-n) = (68) 

We can relate statistical expectation values of the process Ut to some expressions 
from an operator theory. Applying cq.(67) we obtain 

/ duexp{—Q{u)){exp{—TH)exp{—Q)(j)){u) = J duexp{—2Q{u))E[(l){ur{u))] 

(69) 

In the Jiittner model 

V = ^7' + lj^f3^c^ exp{2u) - 1(1 + 37')7'/3cexp(u) (70) 
The eigenstate of H with the zero eigenvalue is 

exp(-n) = exp {^u- ^/?cexp(u)) (71) 
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The statistical weight in eq.(69) 



exp{-2Q{u)) = exp(-/3cr) = 



EI 



(72) 



is just the Jiittner distribution. It can be seen from eq.(69) that in general 



for any equilibrium distribution ^ei- 

The potential V has a minimum. There is a well which supports some bound 
states of H. Eq.(72) defines one of these (normalizable) bound states. All the 
eigenstates of the potential (70) are known ( this is the soluble Morse potential 
[40]) . On the set of functions 



{n a natural number) the expectation value (67) could be calculated explicitly 
using the solution (60). The formula 



allows to calculate the eigenvalues ek and eigenfunctions 0fc from eqs.(60) and 
(67). Wc obtain in this way tlic^ well-known cngenfunctions and eigenvalues of 
the Morse potential [40]. However, for an arbitrary friction we would apply the 
relation (63) the other way round exploiting quantum mechanical methods for 
estimates of the diffusion equilibration. 

7 Discussion:The problem of time 

In general, we have two candidates for time: r and In the relativistic 
dynamics, from the zero component of eq.(ll), we can express r by Hence, 
we can consider the evolution as a function of r or as a function of Let 
us discuss this problem for the diffusion equation (22) assuming (for simplicity 
of the argument) that depends only on the momenta. First, we find the 
solution (x(r, x, p), p(t, p)) of the stochastic equations 



exp(-2n) = 




n 




(73) 



k 



dpj = + ^j'^PjPkd^ In^ERdr + '-fPjedh{T) 



(74) 



dx^ = Pjdr 
dx° = -\p\dT 



(75) 
(76) 



with the initial condition (a;,p). Then, 




(77) 
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is the solution of eq.(22). 

We can obtain from For this purpose integrate both sides of eq.(15) 
in the interval [0,T]. We have 



-1-1 



HniT^ooT-' / dr$, = <S^r (78) 
Jo 

because the Ihs of eq.(78) satisfies the equation Q*^r = 0. Hence, if we denote 

*r = ^ER^T (79) 
Then, the r-average of satisfies the transport equation (16) 

\p\do^R = ^j^PjPkd^d'^^R + 2j^p,d^^R + ^i^PjPkid^ ln$Efl)a'=*fl 

(80) 

We can obtain a solution of cq. (80) from cq. (22) by a random change of time. 
Let us treat the formula (49) for x° (r) as a definition of r 



Jo 



ds (81) 

(here Ps is the solution of eq.(74); x'^{t) can be defined implicitly by (81)). We 
can see from eq.(81) that r depends only on events earlier than As a con- 
sequence pj{x^) = pj(x'^(r)) is again a Markov process . Then, differentiating 
the momenta and coordinates according to the rules of the Ito calculus [28] we 
obtain the following Langevin equations (for mathematical details of a random 
change of time see [39] [28]) 

dpj{x'>) = 2j^\p{x")\-^pjdx^ + l^^\p{x°)\-^pjpkd^ In^ERdx^ 
+'ypj\p{x°)\-iedh{x°) ^ ' 

dx^ =Pj\p{x^)\-^dx^ (83) 

Let \l/(x, p) be an arbitrary function of x and p and (a;(x", x, p),p(a;°,p)) the 
solution of eqs.(82)-(83) with the initial condition (x, p) then 

(a;0 , X, p) = [vt (x(a;0 , X, p) , p(a;0 , X, p) ) ] (84) 

is the solution of the transport equation (80) with the initial condition ^'(x, p). 
We can prove that cq.(80) is satisfied by differentiation of eq.(84) applying 
the rules of the Ito calculus (or the well-known relation between Langevin 
equation and the diffusion equation). Hence, in principle, having a solution 
E[^(xT,Pr)] of the diffusion equation (22) as a function of r we can obtain a 
solution E[^{x{x°),p{x°))] of the transport equation (80) as a function of x°. 
This is like in the deterministic case. The difference is in the averaging over 
events expressed by the expectation value (84) ( a choice of the time parameter 
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and a random change of time in a relativistic diffusion is also discussed in [41]). 
Surprisingly, the new diffusion obtained in this way is the same as the one de- 
rived by the averaging over t. If the limit t ^ cxd of exists then this limit 
is equal to the one resulting from the averaging (78) (i.e., ^"00 = '^r)- We can 
show this property using the equivalence to quantum mechanics of sec. 6 . We 
apply the expansion (73) of ^j- and perform the r averaging term by term. We 
can see that 



exists and is different from zero if and only if Cfe = 0. Hence, the averaging over 
T gives the same result as the limit r — » 00 of ^fr- In this sense for large r the 
probability distribution ^'t- depends only on x^. 

If the initial distribution $ depends on x° then we must perform the limit 
T — > 00 (or averaging) in order to get rid off the additional time in However, 
if $ docs not depend on .t" (what is a reasonable choice of physical initial phase 
space distribution) then apart from the procedure discussed above which leads 
from (x, p) to $ij(x'', X, p) there is still the question whether $t(x, p) could 
possibly have a physical meaning (the solution of the Kompancetz equation is 
of this type). It satisfies eq.(22) instead of the transport equation (80). The 
difference in the r and evolutions is in the Pq^ factor in the collision term 
of the Boltzmann equation. In the massive case po — mc in the non-relativistic 
approximation. Hence, the difference between r and may be insignificant. 
However, in the massless case po — |p| can never be approximated by a constant. 
For a small time r ~ |p[~^x°. We can make this replacement in all solutions 
of proper time equations (in particular in the Kompaneetz equation). However, 
in general, we should make the random time change discussed in this section. 
Nevertheless, for an estimation of the Sunyaev-Zeldovich CMBR spectrum dis- 
tortion [33] [36] [34] this problem is not relevant because the time r is small and 
chosen as a physical parameter (plasma penetration depth) . In such a case the 
spectral deformation is ns + TQ*nE (as discussed at the end of sec. 5) 

It can be seen from our discussion that only a small part of the theory 
of relativistic Brownian motion is applied in astrophysics. The behavior of the 
diffusion at large time (apart from its final effect: the equilibrium) seems to have 
no observational consequences. However, with the increasing sophistication of 
the astrophysical observations the large time effects of Compton scattering may 
give important information about the sources of the the CMBR distortions. In 
this paper we have completely neglected the effect of the gravity. The relativistic 
diffusion of massive particles on a general gravitational background is discussed 
in [42] [43] [44] [45] . We elaborate relativistic diffusion with a friction leading to 
equilibration in a gravitational field for massive as well as massless particles 
in a forthcoming paper [46]. In such a framework the role of strong gravity 
(resulting from compact stars or dark matter) and the effect of the expansion 
of the Universe upon the photon diffusion can be investigated. 



lim T 



1 
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